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Abstract. Using a new long X-ray observation of the cluster of galaxies Sersic 159-03 with XMM-Newton, we derive radial 
temperature and abundance profiles using single- and multi-temperature models. The fits to the EPIC and RGS spectra prefer 
multi-temperature models especially in the core. The radial profiles of oxygen and iron measured with EPIC/RGS and the line 
profiles in RGS suggest that there is a dip in the 0/Fe ratio in the centre of the cluster compared to its immediate surroundings. 
A possible explanation for the large scale metallicity distribution is that SNIa and SNII products are released in the ICM through 
ram-pressure stripping of in-falling galaxies. This causes a peaked metallicity distribution. In addition, SNIa in the central cD 
galaxy enrich mainly the centre of the cluster with iron. This excess of SNIa products is consistent with the low 0/Fe ratio we 
detect in the centre of the cluster We fit the abundances we obtain with yields from SNIa, SNII and Population-Ill stars to derive 
the clusters chemical evolution. We find that the measured abundance pattern does not require a Population-Ill star contribution. 
The relative contribution of the number of SNIa with respect to the total number of SNe which enrich the ICM is about 25- 
50%. Furthermore, we discuss the possible presence of a non-thermal component in the EPIC spectra. A potential source of 
this non-thermal emission can be inverse-Compton scattering between Cosmic Microwave Background (CMB) photons and 
relativistic electrons, which are accelerated in bow shocks associated with ram-pressure stripping of in-falling galaxies. 
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' 1. Introduction 



Hot diffuse X-ray emitting gas dominates the visible mass 
in clusters of galaxies, but the structure and evolution of the 
cluster is not yet fully understood. During its formation, su- 
pernova explosions and galactic winds of member galaxies 
have enriched the Intra-Cluster Medium (ICM) substantially 
(De Young 1978). The abundances can provide the relative 
contribution of Supernova type la (SNIa), Supernova type II 
(SNII) and population-III stars (PopIII) to the enrichment of 
the ICM (e.g Iwamoto et al. 1999; Tsuiimoto et al. 1995), be- 
cause the abundance ratios of the various elements are signa- 
tures of supernova types la and II and possibly of the remains 
of PopIII stars (e.g. Gibson et al. 1997; Loewenstein 2001; 
FBaumgartner et alJl2005l) . The radial distribution of the metals 
provides information about dynamical ways to enrich the ICM 
like, for example, ram-pressure stripping (e.g. Schindler et al. 
l2005l) . 

Clusters of galaxies appear as knots in the cosmic web. 
They accrete gas from the surrounding filaments, which con- 
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sist of warm gas with a temperature in the range of 10^ ^ K. 
According to numerical hydro-dynamical simulations by e.g. 
Cen & Ostriker ( 1999) and Dave et al. (2001), this Warm-Hot 
Intergalactic Medium (WHIM) could contain about half of the 
missing baryons in the universe, 

There have been sever al attempts to det ect th is WHIM in 
emission. In the late 90's iLieu et all ^§9^ and iMittaz etall 
( 1998) discovered a soft X-ray excess in EUV and ROSAT 
spectra of several clusters. More recent observations with 
XMM-Newton appear t o confirm the presence of a soft excess 
in some cluster spectra ("Kaastra et al."2003^. 'pinogue nov et all 
2003). The detection of a possibly redshifted O VII line, which 
traces gas with a temperature of ~ 10* K, suggests that the gas 
might be the WHIM. Unfortunately, current instruments do not 
have sufficient spectral resolution to prove that the emission is 
indeed extragalactic. 

Recent measurements of Sersic 159-03 bv lBonamente et all 
(2 005) and Kaastra et al. (2003) show that the soft-excess can 
be fit both using thermal and non-thermal models. Inverse- 
Compton scattering of CMB photons with relativistic elec- 
trons can also contribute a non-thermal power-law com- 
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ponent to the spectra . Th i s me chanism was akeady pro- 
posed by Isarazin & Lied ( Il998l) to explain the extreme- 
ultraviolet emission from clusters. In the hard X-ray band 
up to 80 keV, detections of non-thermal emission have been 
claim ed in several clusters, for exam ple Coma and Abell 
2256 dFusco-Femiano et alJl999ll200 5'). But these BeppoSAX 
detections of the hard-e xcess are still subject to debate 
jRossetti & Molendil2004l) . 

The cluster of galaxies Sersic 159-03, also known as ACO S 
1 101, was discovered bv iSersic. (.1974.) . Since then it was stud- 
ied in X-rays as part of several cluster samples: e.g. EXOSAT 
(Edge & Stewart 1991) and ROSAT dAllen & Fabian 1997). 
iKaastra et al. i2001) reported results from an XMM-Newton 
observation with a useful exposure time of about 35 ks show- 
ing a radial temperature profile which peaks at kT - 2.7 keV at 
a radius of 2' from the core. The temperature drop in the core 
is relatively modest, while the temperature outside the 2' radius 
drops rapidly to values around 0.5 keV. Because Sersic 159-03 
is thought to show a large so ft X-rav exce s s, the cluster was 
also included in the sample of lKaastra et alJ ll2003l) . 

In this paper we present results from a 121 ks long XMM- 
Newton Uansen et al.„2001.) observation of Sersic 159-03. The 
main goal of the paper is to obtain accurate temperature and 
abundance profiles as far out from the core as possible us- 
ing a new method for handling the background. We exploit 
the large effective area of XM M to obtain tempe rature and 
abundance profiles wit h EPIC ( Turner et al. 2001) and RGS 
dden Herder et aljf2001» . This deep observation also allows us 
to study the nature of the previously detected soft excess in 
more detail, and in addition it provides more accurate radial 
profiles of the temperature and metal abundances. We then fit 
the derived abundances to yields of Supernovae type la and II, 
and PopIII stars. Moreover, we discuss the potential presence 
of non-thermal emission in the cluster. 

Throughout this paper we use Hq = 70 km s ' Mpc ' , On, = 
0.3 and Qa = 0.7. Using this cosmology 1' is 73 kpc at the clus- 
ter redshift of 0.0564 (.Maia et al 1987,,) . The elemental abun- 
dances presente d in this paper a re given relative to the solar 
abundances from lLodder sl (l2003h . 

2. Observations and data analysis 

The XMM-Newton observation of Sersic 159-03 was per- 
formed on November 20, 2002 and had a total duration of 121 
ks. The two EPIC MOS cameras were operated in Full Frame 
mode and the EPIC pn camera in Extended Full Frame mode. 
For all EPIC cameras the thin filter was used. The RGS instru- 
ments were operated in the standard spectroscopy mode. All 
data were analysed with the 6. 1 .0 version of the XMM Science 
Analysis System (SAS). 

2.1. EPIC Analysis 

One of the most important things to account for in extended 
source analysis is the background. Because we intend to mea- 
sure the cluster properties as far from the core of the cluster 
as possible, we need an accurate estimate of the local back- 
ground, especially for the dim outer parts of the cluster. In gen- 



eral the source fills the entire field of view, which makes a direct 
measurement of the local background very difficult. Common 
practice is to extract spectra from a combined event list of sev- 
eral observations of empty fields, like the ones compiled by 
Lumb et al.l (l2002h or Read & Ponman (2003), and use them as 
best estimate for the local background. The spectra extracted 
from this blank field are scaled and then subtracted from the 
source spectra. This method works fine in areas in which the 
surface brightness of the source is high, and with clusters that 
have similar background conditions (e.g. particle background, 
local cosmic background and instrumental background) com- 
pared to the blank fields. 

The Cosmic X-ray Background (CXB), however, con- 
sists of multiple components, which are more extensively de- 
scribed in Sect. 12.1.31 The emission can be roughly divided 
in two parts: a soft thermal component originating from hot 
plasma in our own galaxy and a non-thermal power-law com- 
ponent caused by unresolved distant point sources, predom- 
inantly AGN. The soft galactic CXB component varies spa- 
tially across the sky. Therefore, the total photon background at 
low energies can yarv u p to 30-60% from pointing to pointing 
dRead & Ponmanlf2003l) . It is clear that when the background 
conditions in the source observation are very different from the 
average blank-sky background, this can lead to systematic un- 
certainties in the fitted parameters from the outer parts of the 
cluster. 

In Fig. [Owe show the systematic effect of background sub- 
traction on a temperature profile extracted from the EPIC data 
of Sersic 159-03. For this plot we deliberately introduced an 
error of -1-10% and -10% in the background normalisation. The 
crosses show a profile derived with no background scaling at 
all, while the circles and triangles show the temperature profile 
for a scaling of -10% and -1-10%, respectively. From this plot 
we see that overestimating the background results in a lower 
temperature, and vice versa. Note that we included a 10% sys- 
tematic error on the background during fitting. We conclude 
that the uncertainties in the temperature are mainly caused by 
the background scaling and therefore we can gain a lot in accu- 
racy by constraining the background normalisations. 

In order to find the best estimate for the local background 
we divide the background into three components: soft-protons, 
instrumental background, and the CXB. By separating these 
components in our analysis, we will be able to constrain the 
normalisation and accuracy of each com ponent. The method 
we u se is similar to the one used by IPe Luca & Molendi 
('2004), but with a few modifications. The steps are described 
below. 

2.1.1. Soft Protons 

In order to minimise the effect of soft protons in the spectra, 
we cut out the time intervals whic h show enhanced soft-pr oton 
flux using the method described in Pratt & Arnaud '(■20021). We 
make a light curve using the events with energies between 2-12 
keV and with a bin size of 100 s. From this light curve we ex- 
tract the distribution of counts per bin. The resulting histogram 
is then fitted with a Poissonian function to determine the aver- 
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Table 1. The closed filter observations that we use for the EPIC 
background subtraction. 
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Fig. 1. Systematic effects in a temperature profile obtained with 
a blank-sky background dataset (Read & Ponman 2003). The 
points with the cross are derived from spectra with no back- 
ground scaling. The background for the circled data points is 
scaled with 0.9 and the triangles with 1.1. 



age number of counts per bin. The thresholds are subsequently 
fixed to ± 2 V^, where is equal to the mean number of 
counts within a 100 s bin. This way we cut out every flare 
which reaches the 2cr level. This method is more strict in re- 
jecting soft-proton flares than using the conventional >10 keV 
light curve with 3cr clipping. Because of the broad energy in- 
terval we choose, the signal-to-noise in the light curve is very 
high and even small flares stand out from the quiescent level. 
By putting a strict 2cr threshold we exclude also small flares 
without losing a lot of usable data (~ 4% max. is lost in this 
case). This way we minimise the effect of soft-proton flares 
which can bias our temperature determination like in Fig. ^ 
After applying the threshold we obtain an effective exposure 
time of 81 ks for both MOS instruments and 60 ks for pn. 

The disadvantage of this method is that we do not necessar- 
ily subtract the quiescent level of soft protons. We can get an in- 
dication of the magnitude of the soft-proton count rate by using 
CCD9 of the RGS instrument. For this CCD the photon count 
rate from the source is very low, so the count rate is dominated 
by soft protons and instrumental noise. For our observation we 
find that the RGS CCD9 count rate is 0.086 cnts/s in the quies- 
cent periods. Compared to all the soft-proton background count 
rates measured during the lifetime of XMM-Newton, this resid- 
ual soft proton background is low ( Gonzalez-Riestra 2004). 
Therefore, the influence on our analysis should be very small. 



2.1.2. Instrumental background 

A good template for the instrumental background can be ob- 
tained from a closed-filter observation. For this analysis we use 
a closed-filter observation which was performed after the cool- 
ing of the EPIC and RGS instruments in November 2002. The 
properties of the observations are fisted in Tabled 



Instrument 


Obs. date 


Obs. ID 


Exposure 


MOS 


2003-04-06 


0150390101 


200 ks 


pn 


2003-08-12 


0160362801 


35 ks 



The instrumental background contains roughly two compo- 
nents: fluorescenc e lines and a power law. A study of the EPIC- 
pn background bv'Katavama et all (jSooT) shows that the vari- 
ation of the fluorescence lines is much lower than the variation 
in the continuum. We use the power-law component to model 
the variable hard-particle background and the intrinsic instru- 
mental noise. This background component is mainly caused by 
hard-particles which are able to reach the detector even when 
the filter wheel is in closed position. We especially choose a 
long closed-filter observation, not taken during the passage of 
the Earth radiation belts, to avoid getting a too large variation 
in the hard -particle flux with respect to our Sersic 159-03 ob- 
servation. 

Because these variations can have a significant effect on fit- 
ted parameters (see Fig.QJ, we would like to be able to modify 
the normalisation of the hard-particle power law by a few per- 
cent without changing the normalisations of the instrumental 
lines. In order to find the scaling factors for this normalisation, 
we can use the events registered outside the field of view (out- 
of-FOV) of the EPIC instruments in both the closed-filter and 
source observations. These non -illuminated parts of the CCD 
chips provide unbiased data on the power-law normalisation, 
also when a bright source is present in the field of view. 

We start by determining the power-law index of the instru- 
mental power law by fitting a closed-filter spectrum extracted 
from the full field of view. We assume that this power-law index 
will be constant over the whole detector, since the observation 
does not suffer from vignetting and PSF effects of the mirror. 
We fit the spectrum with just a spectral redistribution file and 
no effective area file. From these fits we obtain a photon in- 
dex (F) of 0.15 for MOS and 0.37 for pn. The MOS value is 
consistent with the value o f ~0.2 obtained in the same way by 
iDe Luca & Molendil (l2004 

We extract the events labelled as out-of-FOV events 
(FLAG==#XMMEA_16) to determine the normaUsation of the 
instrumental background. To be sure that photons and soft- 
protons that scatter into the shielded out-of-FOV region of 
the detector are not polluting the measurement, we only count 
events registered outside a radius of 15.4' from the centre of 
the field-of-view. We divide the out-of-FOV count rates in our 
Sersic 159-03 observation with the closed-filter out-of-FOV 
count rates in the 8-12 keV band to obtain an instrumental 
background scaling factor (c). The c values are 1.03 + 0.02, 
0.97 ± 0.02 and 0.97 + 0.06 for MOS I, MOS2 and pn respec- 
tively. In this case the correction on the power-law normalisa- 
tion is only ±3% and within the statistical error, but we prefer 
to use the most likely value for c. Despite the fact that our val- 
ues for c are consistent with being 1 , this is a useful exercise, 
because the normafisation of closed filter spectra is known to 
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vary by about 15% JPe Luca & Molendill2004 . Already a de- 
viation of a few percent can lead to discrepancies in the tem- 
perature determination. 

We cannot simply multiply the closed-filter spectra with 
our values for c, because then the instrumental lines would not 
be well subtracted. The only component that we have to scale 
is the power-law component. Therefore, we add (or subtract) a 
power law with the same slope (F) to the original closed filter 
spectra (S'oiig). The normalisation of the power law (n) we add 
is derived as follows. We first solve this simple system forS'oiig 
to approximate the scaled spectrum 5 seal : 

seal — orig ^ j 

S seal — nE + S orig , 

and work out n; 

— orig ■ 

We can subsequently calculate n for every energy bin of energy 
E. We exclude energy intervals affected by instrumental lines. 
The mean value of n over all these bins is determined by fitting 
a Gaussian to a histogram of n. This value will be used to add 
or subtract the power law from the closed filter spectra. Finally, 
we subtract this scaled spectrum from the source spectra. 



Table 2. The CXB background components we use in the fits of 
Sersic 159-03. The integrated absorbed intensity over the 0.3- 
10 keV energy range was calculated using an A^h value of 1.79 
X 10^'^ cm"^. In the thermal components the abundances are set 
to solar values. 



Component 


kT 


r 


z© 


Integrated intensity 




(keV) 






(erg cm"- s"' deg"^) 


Soft thermal 


0.070 




0.3 


2.23 X 10-'^ 


Hard thermal 


0.20 




1.0 


1.03 X 10^" 


Power law 




1.41 




3.16 X 10 " 



Table 3. Boundaries of the annular extraction regions used in 
our EPIC analysis. The annuli are centred on the cluster centre. 



Annulus 


Inner boundary (') 


Outer boundary (') 


1 


0.0 


0.5 


2 


0.5 


1.0 


3 


1.0 


2.0 


4 


2.0 


3.0 


5 


3.0 


4.0 


6 


4.0 


6.0 


7 


6.0 


9.0 



2.1 .3. Cosmic X-ray background 

Instead of subtracting the CXB from the spectra, like it is done 
with blank-field data, we include the CXB components during 
fitting. Because the CXB can vary spatially across the sky we 
can more easily adapt the flux of each component to the local 
conditions. Unlike the case of the instrumental background and 
soft-protons, all the photons enter through the XMM mirrors in 
the same way as the source photons. Therefore, the response 
files (rmf) and effective area files (arf) must be applied and we 
must fit the CXB simultaneously with the source spectra. 

For the fitting we use the s pectral components describe d by 
iKuntz & S nowden'('2000') and be Luca & M olendi'('2004'). For 
the thermal components in this pa per we use MEKAL mod- 
els, but the thermal components in iKuntz ^ ^ nc^wc^^ij 
are fi tted with Raymond-Smith models ( Ravmond & Smith' 
Il977l) . There are, however, significant differences in ionisation 
balance and line strengths between MEKAL and Raymond- 
Smith, because the MEKAL code includes many more ions 
and up-to-date atomic constants. Hence, a temperature deter- 
mined from a spectrum using the Raymond-Smith model is 
different from a temperature determined using MEKAL for 
the same dataset. To obtain consi stency, we choose t o con- 
vert the model temperatures from iKuntz & SnowdeJ ll200(]|) 
to MEKAL temperatures to match the measured spectra by 
iKuntz & Snowd^n ( 2000) in our model ling. Unfortunately, this 
transformed Kuntz & SnowdenI (l2000l) model does n ot fit a 
spectrum extracted from the 'Read & Ponman' ("2003") blank- 
fields. Therefore, we empirically fit the Read & Ponman ( 2003) 
blank-fields with a power law, a soft, and a hard thermal com- 
ponent, which turns out to be a good description of the data. 
Based on the best fit and the model by .Kuntz & Snowdeni 



(1200(1) we fix the temperatures and metallicities to the values 
listed in Table|2 

In order to estimate the normalisations of the local back- 
ground components around Sersic 159-03 we fit them to the 
outermost annulus (9-12' from the core of the cluster) u sing the 
temperatures derived from the iRead & PonmanI (120031) blank- 
fields. We add an additional thermal component for the cluster 
emission. In our final fits for all annuli we fix the normalisa- 
tions of the background components to the fitted values of the 
9-12' annulus, which are listed in Table|2l For comparison, the 
2-10 keV absorbed integrated intensity of the power-law com- 
ponent is 2.26 X 10"" erg cm"^ s"' deg"^, which is consistent 
with the value of (2.24 + 0.16) x 1 0~" erg cm"^ s ' deg"^ 
found bv lDe Luca & Molendil(l2004 . 

In order to avoid large background fluctuations in our ex- 
traction regions we cut out bright point sources which were 
identified by eye. The unresolved point sources are taken into 
account by the CXB power law. 

2.1 .4. Data extraction and fitting 

As we are interested in the spectral parameters as a function of 
radius, we extract spectra from annuU as defined in Table|3l To 
gain statistics we choose a width for the annuli which is > 30". 
By using these relatively wide annuli, we are also less sensi- 
tive to the energy dependent shape of the PSF and therefore we 
neglect this effect in the rest of the analysis. Because of small 
calibration differences between MOS and pn we include a 5% 
systematic error on the source and background spectra. This 
error is added in quadrature to the poissonian error of the data 
points. 
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Fig. 2. RGS extraction regions projected on top of a MOSl image (top panel) and projected on top of a combined RGS1/RGS2 
detector image (bottom panel). In the eighth and ninth CCD of RGS the instrumental effect of fixed-pattern noise causes the 
fish-bone shaped patterns. 



2.2. RGS analysis 

W e extract the R GS spectra following the method described 
in iTamura et alJ (^001). In addition, we do spatially resolved 
spectroscopy using the RGS data: with RGS it is possible to ex- 
tract spectra from different regions in the inner 4' of the cluster. 
The different regions can only be separated if they lie along the 
cross-dispersion axis of the instrument. We select the events 
from several rectangular areas on the CCD strip in the cross- 
dispersion direction. In Fig. |2we show how these extraction 
regions are projected on the sky (MOSl, upper panel) and on 
the chips of the RGS (lower panel). 

Because the cluster fills the entire field-of-view of the RGS, 
we need a blank field observation to extract the background 
spectrum. For this observation, which was taken prior to the 
cooling of RGSl but after cooling of RGS2, we choose a 
Lockman Hole observation with an effective exposure time of 
100 ks performed just 3 XMM-Newton orbits after our obser- 
vation, also still before the cooling of RGS 1 . This way we min- 
imise systematic effects due to hot pixels and instrument re- 
sponse. The flare subtraction is analogous to the method used 
with EPIC, but now we use the events from CCD 9 outside the 
central area. We select only events that have a position greater 
than 30" from the dispersion axis to make the light curve. This 
method was applied to both source and background datasets. 

The RGS spectrometer operates without a slit. This means 
that all photons from within the (in our case) 5' x ~ 12' field 
of view end up in the final spectrum, but not necessarily at the 
right wavelength. Only photons which are emitted in the cluster 
centre end up at the dispersion coordinate which corresponds 
to the coiTect wavelength. If the photon originates from the out- 
skirts at an angle 8 (projected on the dispersion axis) from the 
cluster centre, then the instrument will register it at a different 
dispersion coordinate and assign a wavelength to it which is 



shifted with respect to the true wavelength. This shift in wave- 
length depends linearly on incidence angle (arcmin) projected 
on the dispersion axis: 



A/l = 0.138AA6i. (3) 



Because of this effect, the line-emission appears to be broad- 
ened depending on the spatial extent of the source along the 
dispersion direction (see .Davis.2001. for a complete discussion 
about grating responses). 

In order to describe the data properly, the spectral fits need 
to account for the effect described above. In practice, this is ac- 
complished by convolving the spectral models with the surface 
brightness profile of the source along the dispersion direction 
(Tamura et al. 2004). Therefore, we extract the cluster inten- 
sity profile from MOSl along the dispersion direction of RGS. 
For each extraction region we can convert the spatial profile, 
which is a function of 6, into a line profile in A using Eq.|3] We 
convolve this profile with the model spectrum during spectral 
fitting. However, the width of the measured line profile can be 
different from the line profile that we derive from MOSl, be- 
cause the surface brightness at the wavelength of a spectral line 
is not necessarily the same as the surface brightness in a broad 
energy band. To take this effect into account, we multiply the 
wavelength axis of the derived line profile with a scale factor A 
scale factor of 1.0 corresponds directly to the wavelength scale 
of the profile we derive from the continuum. A scale factor of 
2.0 stretches the line profile in wavelength space by a factor of 
two compared to the original line profile. In this way, we are 
able to fit the widths of the line profiles and to have a measure 
of the spatial extent of the line-emission region. 
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3. Spectral Models 

In our analysis we fit several models to the spectra using 
the SPEX package (Kaastra et al. 1996). These models can 
be a combination of a number of thermal models (MEKAL) 
and a power-law model. Two models, however, include a 
more sophisticated combination of thermal models that we call 
Diffe rential Emission Measure (DEM) models. P r evious papers 
(e^lPeterson eT7\ Eooi iKaastra et al.l Eooi Ide Plaa et al.l 
I2OO4I) show that many clusters can be fitted better when we 
use a distribution of temperatures instead of a single temper- 
ature model. Actually, we expect that the plasma within our 
annuli should contain multiple tempe ratures and not one. An 
observational proof of this is given in lWerner et alJ ll2006h . In 
the case of Sersic 159-03 we fit two types of distributions: a 
truncated power-law distribution (wdem) and a Gaussian distri- 
bution igdem). 

3.1. WDEM 

We use the so-called wdem model, where the emission mea- 
sure, Y - J rienndV, of a number of thermal components is 
distributed as a truncated power l aw. This is shown in Eq. @ 
adapted from lKaastra et alJ ll2004 : 

dY _ ( cT'/" /3T^,, <T < 

df~\0 T>T^,,WT </3T^,,. ^ ' 

This distribution is cut off at a fraction of Tmax which is /3T^^^. 
The value of fi is set to 0. 1 in this study. The model above is an 
empirical parametrisation of the DEM distribution found in the 
core of many clusters. In this form the limit ff — » yields the 
isothermal model at Tmax- 

In order to compare the outcome of the wdem model 
with single-temperature models, we can calculate the emission 
weighted mean of the DEM distribution. The mean temperature 
^^niean follows from Eq.|5] 

fT§dT 



• mean r 

J§dT 

When we integrate this equation between j3T„ 



(5) 



and Tmax, we 



a and /?: 



(6) 



we present in this paper are calculated 



(2-Hl/a')(l 
The values for fcrmean 
using Eq.|6l 

A detailed comparison of the wde m model with t he cla ssi- 
cal cooling-flow model can be found in lde Plaa et all tOOA In 
general the wdem model contains less cool gas than the classi- 
cal cool ing-flow model, which is c onsistent with recent obser- 
vations JPeterson et alJ200lLEo03l) . 



In this equation x - log T and xq - log Tq where Tq is the aver- 
age temperature of the distribution. The width of the Gaussian 
is cTj. Compared to the wdem model this distribution contains 
more emission measure at higher temperatures. 

4. Results 

4.1. EPIC/RGS spectral fitting 

In this section we fit both the EPIC and RGS spectra with the 
models described in Sect. |3] From the fits we obtain temper- 
atures and abundances for which we also derive spatial infor- 
mation. In all the fits we fix A^h t o 1.79 x 10^" cm~^ which 
is the value deduced from Hi data dDickev & Lockmanl l99flt). 
This value is the same as the one used in*Ka astaetaLn2004B . 
We use the Verner et al. (1996) cross-sections in our absorp- 
tion model. Throughout the paper we use errors at the A^^ = 1 
(68%) level for one interesting parameter 

The EPIC spectra are rebinned to the optimal binning of 
~l/3 FWHM of the resolution of the instrument. We use the 
response matrix to derive the FWHM for every energy. The 
spectra are fitted over the 0.3-10 keV range. We are aware that 
below 0.6 keV there are some calibration uncertainties of the 
order of 5-10% when using SAS 6.1.0. This is partially ac- 
counted for by the 5% systematic error that we add to the spec- 
tra during fitting. The remaining caUbration uncertainties are 
still small compared to the uncertainties in the soft X-ray back- 
ground. Because of the importance of the soft-excess problem 
and the role of O VII we choose to include the 0.3-0.6 keV 
energy band in our fits. 

For RGS we discuss the spectrum which was extracted 
from a 4' wide strip in the cross-dispersion direction, with 
which we obtain the highest signal-to-noise. In addition we also 
present spatially-resolved spectra extracted from smaller strips 
in the cross-dispersion direction. This provides a radial profile 
of the temperature and abundances in the core region (0-2'). 

After preliminary analysis, the RGS spectra show a discrep- 
ancy with any thermal model around 29-33 A. The data points 
are significantly below the models in this interval, as can be 
seen in Fig. |5] Since the wavelength range coincides with a 
read-out node of CCD2, we believe this feature is instrumental 
in nature. If we perform the same analysis on an earlier ob- 
servation of Sersic 15 9-03 the discrepancy disappears, which 
supports our view. In IPeterson et al.l (I2OO3I) some other RGS 
observations of clusters of galaxies show the same problem, 
namely Hydra A, Abell 496 and MKW 3s. A large study of 
many RGS observations shows that about ~ 5% of the observa- 
tions show an anomaly in the count rate of CCD2 (A. Pollock, 
priv. comm.). We therefore believe the feature we observe is 
instrumental in nature and we ignore the 29-33 A interval in 
our analysis. 



3.2. GDEM 

Another DEM model that we use is a Gaussian differential 
emission measure distribution, gdem, in log T: 

Y(x) = ^£_e-(-'^--*-o)'/24. (7) 



4.1 .1 . Integrated spectrum within 4' 

In Table|4]we show the results from fits to a spectrum extracted 
from a circular region which has a radius of 4' and is centred on 
the core. In this circular region the cluster signal is well above 
the background. We can use the high signal-to-noise ratio of 
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Fig. 3. EPIC spectrum extracted from a circular region of 4' around the core. The total best-fit model, the wdem model component 
and the background component are shown. 



this spectrum to obtain accurate values for all the fitted parame- 
ters. The;(f^ value of the single-temperature fit (1228/916) indi- 
cates that a single-temperature model is potentially not the best 
description of the data. Multi-temperature models, however, 
do produce a. x^- close to 1.0. Hence, our abundance analysis 
will be based on the multi-temperature models. There are other 
models that fit the data equally well. These alternative models 
will be discussed in Sect. 14.31 The potential presence of a rel- 
atively small contribution of other emission components how- 
ever is of no consequence for the observed trends we present in 
this section. 

An example of a spectrum fitted with a DEM model is 
shown in Fig. |5] When we set the line emission to zero in our 
model and plot the residuals, we obtain the plot shown in Fig.|3] 
The lines from all elements for which we fit abundances are ev- 
ident in the spectrum. Especially the iron lines are very strong, 
but also silicon, sulfur and calcium have a high signal-to-noise. 
Oxygen, magnesium, argon and nickel are clearly visible: this 



plot shows the total spectrum; the strengths of these lines are 
weaker in individual annuli. 

The spectral models constrain the abundances very well, 
except for neon and magnesium, because these lines are 
blended with iron at this spectral resolution (see Fig. 0}. The 
determination of the oxygen abundance becomes difficult in the 
outer parts of the cluster Here, the oxygen in the galactic fore- 
ground emission starts to play a more important role and could 
bias the abundance measurement. In the central region (< 4') 
the cluster flux is high enough to get an accurate measurement. 
In general, all abundances are between 0. 1 and 0.6 times their 
solar value (Table 0}. Only neon shows high values up to 1, 
but since the lines of this element are blended with the Fe-L 
complex, the value is highly correlated with the temperature 
distribution. 

In Fig.|5]we present the total RGS spectrum of Sersic 159- 
03 extracted from a 4' wide strip in the cross-dispersion direc- 
tion and centred on the core. Because of the high statistics, we 
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Residuals of a wdem fit to the EPIC total spectrum extracted from a circle with a radius of 4'. The line emission in 
is set to zero to show the line emission on top of the continuum. The left panel shows the residuals from 0.3-5 keV 
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Table 4. Fit results for an EPIC spectrum extracted from a cir- 
cle with a radius of 4' and centred on the core. Fluxes are cal- 
culated over the 0.3-10 keV range and presented in 10"'° erg 
cm"^ s"' deg"^. Emission measure (^thermal = JnenndV) is 
given in 10^^ cm"-' and Fpow is given in 10^' ph s"' keV"' at 1 
keV. 



Parameter 


single-temp 


wdem-model 


gdem 


■model 


^thermal 


20.31 ±0.15 


20.55 ±0.15 


21.07 


±0.15 


^thermal 


19.26 ±0.14 


19.40 ±0.14 


19.46 


±0.15 


kT 


2.568 ± 0.009 




2.472 


+ 0.010 


kTiYisun 




2.60 ± 0.03 






kT 

"■^ max 




3.41 ± 0.03 






a 




0.45 ± 0.02 






(Tt 






0.226 


± 0.005 





0.36 ± 0.05 


0.30 ± 0.04 


0.19 


± 0.03 


Ne 


0.89 ± 0.08 


0.16 ± 0.08 


0.000 


± 0.013 


Mg 


0.02±0.04 


0.11 ± 0.04 


0.08 


±0.04 


Si 


0.196 ±0.017 


0.239 ± 0.017 


0.252 


±0.018 


S 


0.122 ± 0.017 


0.18 ±0.02 


0.20 


± 0.02 


Ar 




0.14 ± 0.04 


0.19 


± 0.05 


Ca 


0.31 ± 0.05 


0.36 ± 0.05 


0.44 


± 0.06 


Fe 


0.360 ± 0.007 


0.346 ± 0.008 


0.242 


± 0.005 


Ni 


0.42 ± 0.09 


0.45 ± 0.08 


0.35 


± 0.08 


x'/dof 


1228/916 


1017/915 


949/915 



Table 5. Fit results for RGS spectra extracted from a 4' wide 
strip and fitted over a 8-38 A range ignoring data from CCD2 
(see text). Emission measures (J' = / n^nndV) are given in 
lO^*" cm"^. The iron abundance is given with respect to solar 
abundances. The scale factors for oxygen and iron mentioned 
below are those explained in Sect. 12.21 



Parameter 


single-temp 


wrfem-model 


grfem-model 


Y 


7.27±0.06 


7.08±0.07 


7.25±0.06 


kT (keV) 


2.6±0.09 




3.25±0.18 






3.04 ±0.16 




kT„,,, (keV) 




4.0±0.2 




a 




0.46±0.03 




a 






0.28±0.02 


C/Fe 


0.00±0.19 


0.0±0.4 


0.0±0.3 


N/Fe 


0.00±0.18 


0.1±0.5 


0.0±0.5 


0/Fe 


0.85±0.10 


0.87±0.10 


0.91±0.13 


Ne/Fe 


1.04±0.15 


0.73±0.12 


0.71±0.14 


Mg/Fe 


0.32±0.11 


0.35±0.11 


0.32±0.12 


Fe 


0.98±0.06 


1.20±0.08 


1.14±0.09 


Scale (O) 


1.1±0.2 


1.6±0.3 


1.5±0.3 


Scale (Fe) 


0.35±0.05 


0.36±0.05 


0.31±0.05 




997 / 856 


925 / 854 


963 / 855 



can also use the second order spectrum. In both orders the Fe-L 
spectral line complex between 10-15 A and the Mg line near 9 
A are well resolved. Above 19 A the first order spectrum shows 
a prominent O viii Ly a line. 

The fit results for RGS are shown in Table |5] The elemen- 
tal abundances determined from the different models are con- 
sistent with each other within 3cr. Unfortunately, the carbon 
and nitrogen abundances are not well constrained. The width 
of the oxygen line is consistent with the width of the line pro- 
file derived from the continuum emission in MOSl, and the 
width of the iron hnes is about 1/3 smaller (Fig.|6}. This sug- 



gests that the oxygen distribution is more extended across the 
cluster centre, while the iron abundance is strongly peaked in 
the centre. Looking at the values we see that the fit to the 
RGS data does not strongly prefer a DEM distribution over a 
single-temperature model. All models are acceptable, cont rary 
to the merging cluster 2A 0335-1-096 described bv .Werner et alJ 
(f2006l) . 

4.1.2. Radial profiles using thermal DEM models 

By fitting the spectra extracted from the annuli for EPIC and 
from strips for RGS, we can make radial profiles of tempera- 
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Fig. 5. 1 " and 2'"^ order RGS spectrum of Sersic 159-03 extracted from a 4' wide strip centred on the core. The continuous line 
represents the fitted wdem model. On the x-axis we show the observed wavelength. 



tures and abundances. Again, we fit the spectra with a single- 
temperature model, a wdem model and a gdem model (Table|6j. 
In Fig. (left panel) we present the temperature profile ob- 
tained for these models. For wdem the maximum temperature 
(^T'max) and mean temperature (^rmean) are shown. The profile 
shows a slight increase of the temperature within a 2' radius. 
Beyond 3' the temperature rapidly drops to about 1.5 keV. The 
single-temperature profile is consistent with the profile derived 
from t he earlier XMM-Newton observation by iKaastra et alJ 
1I2OOII). 



The DEM parameters a (wdem. Fig. middle panel) and 
(Tt (gdem, Fig.Qright panel) are quite constant, but increase 
slightly in the outer parts, where the temperature gradient is 
large. This means that the fit needs a broader range of temper- 
atures to fit the spectrum. 

From the fits we obtain radial profiles of the abundances 
of several elements. The radial profiles of oxygen, silicon, and 
iron are shown in Fig.|8l We detect a clear decrease of the sili- 
con and iron abundance with radius. For oxygen the situation is 
less clear: the 0/Fe ratio in the 0-0.5' bin is lower than the ratio 
in the 0.5-4' interval at 3cr confidence level. There is a hint of 



a decrease of the oxygen abundance like we see for iron, but a 
flat distribution cannot be excluded. 

While silicon is quite well constrained independent of the 
DEM model used, the iron abundance derived from the gdem 
model is significantly different from the other two models. The 
absolute values of the gdem oxygen and iron abundances are 
consistent with the ones from silicon. The iron abundance from 
the single-temperature and wdem, however, show a steeper gra- 
dient. From these plots it is clear that the used temperature 
distribution is very important when measuring elemental abun- 
dances. 

An overview of all the fitted parameters is given in Table|6l 
Apart from the elements we discussed in the previous para- 
graph, only sulfur is reasonably well constrained up to 9'. In 
the outermost spatial bin, however, the abundance increases to 
unphysical values. In the core region sulfur shows the same 
profile as silicon. The other elements are less well constrained. 
Neon is always difficult to measure at CCD resolution because 
it is blended with the Fe-L complex near 1 keV. Calcium and 
nickel are only constrained in the core, but the derived values 
are not affected by systematics due to different DEM models. 
In every radial bin the values are consistent with each other. 
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Fig. 7. DEM model fit results for the three models: single-temperature (O), wdem (a) and gdem (*). In the left panel we show 
the radial temperature profiles for all models, including the kTmax (a) and kTme^a (♦) from the wdem parameters. In the middle 
and right panel we show the radial profiles of the DEM parameters a and crj. 




Fig. 8. 
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Abundance profiles of oxygen, silicon and iron derived using three models: single-temperature (O), wdem (a) and : 



The magnesium and argon abundances are poorly determined. 
Only in the core region some points are measured at a signifi- 
cance larger than 2cr. In general these two elements have upper 
limits of about 0.2 times the solar abundances. 

The values for the three models are also shown in 
Table |6] In the core region all fits are acceptable when we use 
the condition < 1.5. However, the xl^ for the single tem- 
peratu re model in the core is much less than in Werner et al.l 
( I2OO6I) in the cluster 2A 0335+096. In the outer parts multi- 
temperature fitting still results in a better fit, but we do not see 
significant differences between the models. 



From RGS we also extract rectangular regions in the cross- 
dispersion direction of the instrument. The regions we use are 
defined in Sect. 12.21 Because the cluster does not show signif- 
icant spatial asymmetries in the cross-dispersion direction, we 
add the spectra extracted from regions with equal distance to 
the dispersion axis: region 1-1-5 and 2 H- 4 as defined in the 
lower panel of Fig. |2l This way we are able to derive radial 
profiles of the core up to a radius of 2'. The best fit values 
for three models are presented in Table The fit results show 
signs of a temperature decrement in the core. However, the tem- 
peratures we determine from RGS are systematically higher 



J. de Plaa et al.: Chemical evolution in Sersic 159-03 observed with XMM-Newton 



11 



Table 6. Fit results for spatially resolved EPIC spectra. The fitted models are single-temperature CIE (s), wdem (w) and a 
Gaussian DEM model (g). Fluxes are calculated over the 0.3-10 keV range and are given in unit erg s ' cm"^ deg"^. Emission 
measures (Y - J n^n[jdV) are given in lO^*" cm"-'. 





M 


0-0.5' 


0.5-1.0' 


1 .0-2.0' 


2.0-3.0' 


3.0-4.0' 


4.0-6.0' 


6.0-9.0' 


Y 




4.75±0.05 


5.55±0.05 


6.17±0.05 


2.88±0.03 


1 .40±0.02 


1 .62 ±0.04 


1.21 ±0.04 






4.72±0.05 


5.58±0.05 


6 93+0 05 


2.91 ±0.03 


1.42±0.02 


1.62±0.03 


1.21 ±0.03 




S 


T-.OOiU.U'-f 






9 97+0 0^ 


1 4zl+n 09 


1 55±0 03 


n 94+0 04 


p 




(2.59±0.02)xlO"" 


(1.038±0.009)xlO"'^ 


(2.76±0.02)xl0"'^ 


(7.72±0.09)xl0""^' 


(3.17±0.05)xl0"'" 


(1.08±0.02)xl0"'" 


(2.59±0.08)xl0"" 




w 


^Z.OU±U.UJ lU 


1 .U4-Z±U.U iU ^X iU 


C9 77-lA (191 V 1 n-9 
\Z.. 1 1 ±U.UZ JX lU 


C7 70_i_n HQAvlO-lO 


ex OA-^A A'ilv 1 A-10 
\J.Z.\j±\j.\jJ ^X 1 u 


Cl CWij-A A1 QAv 1 A-'O 

1 1 .uvj±u.ui y jx lu 


CO 'IQ-lA ASivin-11 




g 


CO -i-H no\N/ 1 


ci n/i"7-i-n nicuvin^s 

1^1 .U4/±U.UiU^X iU 


(Z. /o±U.UZ^X iU 


(, / .£S4±U.Uo}X lU 


(1 0/1 -1- A A'llv 1 A- 10 


C 1 11 0-1- A Al 1 A-iO 
1^1 . i lZ±U.Ulo}X lU 


CO ^Q-lA AOIv 1 A— 1 1 

^z.5y±u.uy^x lu 


kT 




9 ^s^+n 019 
z.jjjiu.uiz 


Z. J J / iv/.v/l J) 


9 c;40+0 01 S 


ftd+Ci 0'^ 


SO+0 04 


1 +0 OS 

Z.Z 1 iV/.U^l 


1 A I -i-A 04 


(keV) 


g 


OQn_i_A All 

Z.ZoU±U.Ul 1 


Z.4/ 1 ±U.Ul J 


O C/1"7iAAlC 

Z.j4/±U.U1j 


CT-lA AO 

Z.J 1 ±u.uz 


/I Q _i_ A A/1 


/1A-i-A AC 
Z.4U±U.U5 


1 Q'3-i-A AA 




w 


1 1A 1 A AT 


T CA 1 A AC 

2.59 ±0.05 


2. 67 ±0.05 


T ^lA 1 A AO 

2.oy±o.06 


2.56±0. 12 


2.27±0. 1 1 




kT-nix 


w 


j.U 1 ±U.Uj 


1 97-^C^ CW 
j.Z / ±U.U'4- 


1 zl9-^0 CK 


J. J i ±U.Uo 


3 31±0 13 


OA-^C^ 10 
Z.VU±U. 1 z 


1 49 A (Mi 
1 .4Z±U.Uo 


a 


w 


0.i5U±0.019 


0.36±0.03 


A ■> A 1 A A"> 

U.39±0.03 


t\ A A I t\ AC 

U.44±0.05 


t\ A'^ I r\ A^ 

0.42±0.07 


0.3 8 ±0.07 


A A 1 1 A AA 

0.01 ±0.09 


CTt 


g 


0.203 ±0.005 


0.208±0.007 


0.209±0.008 


0.236±0.012 


0.244±0.019 


0.24±0.02 


0.14±0.04 


O/Fe 


s 


U.76±U.12 


1.12±0. 13 


1.17±0. 17 


1.3±0.3 


1.2±0.4 


0.2±1.0 


0.0±1.7 




w 


0.76±0.1 1 


1.05±0.14 


1.06±0. 17 


1.2±0.3 


1.1 ±0.4 


0.4±0.6 


0.0±1.7 




g 


0.77±0.1{) 


1.05±0.13 


1.05±0.17 


1.0±0.3 


0.4±0.4 


0.00±0.05 


0.00±0.06 


Ne/Fe 


s 


2.6±0.2 


2.0±0.2 


1.4±0.3 


1.4±0.4 


1.5±0.6 


1.7±1. 1 


0.0±0.9 




w 


0.9±0.2 


0.4±0.3 


0.00±0.15 


0.0±0.2 


0.0±0.5 


0.2 ±0.9 


0.0±0.9 




g 


0.00±0.16 


0.00±0.07 


0.00±0.05 


0.00±0.09 


0.0 ±0.2 


0.0 ±0.5 


0.0 ±0.7 


Mg/Fe 


s 


0.10±0.10 


A AA 1 /\ 1 

0.00±0. 12 


A AA 1 A AO 

0.00±0.08 


0.01 ±0.19 


A l\t\ 1 A AA 

0.00±0.09 


0.1 ±0.5 


0.0±0.4 




w 


0.29±0.10 


0.18±0.13 


0.13±0.15 


0.3±0.3 


0.00±0. 1 1 


0.3±0.5 


0.0±0.4 




g 


0.34±0.14 


0.20±0. 1 8 


0.1 ±0.2 


0.5±0.4 


0.0±0.2 


1.3±0.6 


0.0±0.7 


Si/Fe 


s 


U.56±U.(J5 


A AC 1 /A 

0.45±0.06 


U.57±0.07 


A O O 1 A T 

U.38±0. 12 


A O 1 A AA 

0.43 ±0.09 


0.8 ±0.3 


0.5 ±0.5 




w 


0.69±0.05 


0.56±0.06 


0.70±0.08 


0.51 ±0.1 3 


0.5 ±0.2 


0.8±0.2 


0.4 ±0.5 




g 


1 .07±0.08 


0.88±0.09 


1.08±0.12 


0.8±0.2 


0.9 ±0.3 


1 .4±0.4 


0.7 ±0.6 


S/Fe 


s 


U.4i±U.U0 


0.34±0.07 


0.19±0.09 


0.22±0. 16 


0.3 ±0.3 


0.97±0. 17 


5.6±1.7 




w 


0.56±0.06 


0.47±0.08 


0.35±0.09 


0.40±0.15 


0.5 ±0.2 


1.0±0.4 


5.6±1.8 




g 


0.90±0.09 


0.76±0. 11 


0.57±0.13 


0.7±0.2 


0.9±0.4 


1.3±0.5 


3.5±1.2 


Ar/Fe 


s 


A AT 1 A 1 "5 

0.0z±0.1 J 


0.20±0. 1 / 


0.4±0.2 


0.0±0.4 


0.6 ±0.7 


1.9±1. 1 


16.8±6.4 




w 


0.19±0.14 


0.39±0. 18 


0.7±0.2 


0.2±0.4 


0.8 ±0.7 


1.8±1.1 


17.0±6.4 




g 


0.4±0.2 


0.7±0.3 


1 .2±0.4 


0.4±0.7 


1.4±1.2 


1.6±1.5 


5.4±3.5 




Ca/Fe 


s 


U. jV±U. 1 o 


1 1±0 2 


1 1 ±0 3 


1 0±0 5 


1 6± 1 


3 6±1 5 


jO.O± i J.O 




w 


0.67±0.18 


1.3±0.2 


1.2±0.3 


1.0±0.6 


1.6±1.0 


3.0±1.6 


36.2±15.9 




g 


1.2±0.3 


2.3±0.4 


2.2±0.5 


1.7±1.1 


2.2±1.8 


0.4±2.3 


0.0±6.5 


Fe 


s 


0.453±0.012 


0.381 ±0.010 


0.294±0.009 


0.239±0.011 


0.209±0.013 


0.14±0.02 


0.049±0.009 




w 


0.458±0.011 


0.371 ±0.011 


0.282±0.009 


0.224±0.012 


0.196±0.017 


0.140±0.016 


0.049±0.009 




g 


0.314±0.007 


0.253±0.007 


0.191±0.007 


0.148±0.008 


0.134±0.012 


0.134±0.012 


0.13±0.02 


Ni/Fe 


s 


1.0±0.2 


0.9±0.2 


0.0±0.3 


0.0±0.2 


0.3 ±0.8 


0.0±0.3 


0.0±2.3 




w 


1.2±0.2 


1.0±0.3 


0.2±0.3 


0.0±0.5 


1.1±1.0 


0.0±0.3 


0.0±2.4 




g 


1.4±0.3 


1.1 ±0.4 


0.0±0.5 


0.0±0.8 


2.8±1.5 


0.3±1.6 


6.0±2.5 


x'/dof 


s 


1213/829 


1084/846 


1240/852 


1191/804 


1270/772 


1113/777 


1045 / 757 




w 


1045 / 828 


969 / 845 


1136/851 


1 109 / 803 


1218/771 


1095 / 776 


1045 / 756 




g 


981 /828 


890/845 


1083/851 


1034/803 


1131/771 


1042 / 776 


1079 / 756 



than those from EPIC. This is not surprising, because the RGS 
spectrum also contains emission from the hot gas just outside 
2' which falls within the rectangular field-of-view. From the 
wdem model fits we see that the value of a increases with ra- 
dius, while the width of the Gaussian DEM distribution shows 
a drop in the 1.0-2.0' bin. As expected from the EPIC pro- 
files, the O/Fe ratio is lower in the core than in the outer parts. 
The Ne/Fe ratio, however, is consistent with being flat within 
2' from the core. Finally, the width of the lines, indicated by 
the scale parameter (see Sect.|^) , increases outside the core. 
There it is consistent with being 1.0. In these models we fit 
the average of all the line widths with respect to the continuum 
surface brightness. 

From the O and Fe abundances measured with EPIC and 
RGS in the inner 2' of the cluster, we detect a jump in the 
O/Fe ratios between the inner 0.5' and the annulus from 0.5- 
2.0'. To illustrate this, we combine the single-temperature EPIC 
and RGS results in the 0.5-2.0' region. From EPIC we obtain 
(0/Fe)g 2 o' ~ 1 . 14 ± 0. 10 which is significantly higher than 
the central (0/Fe)gg' = 0.76 + 0.12. The combined O/Fe 
ratio from the 0.5-2.0' RGS results is also significantly higher, 
(0/Fe)gg' 2o' ~ ^-^^ - 0-09, compared to the central value 



(0/Fe)QQ/ Q = 0.53 + 0.05. This jump has a confidence of 
2.5cr and 3cr in EPIC and RGS respectively, which corresponds 
to a combined significance of 3.9cr. Despite the fact that the 
absolute values of the O/Fe ratio are different for the EPIC and 
RGS results, the relative increase of the O/Fe ratio is the same 
for both instruments: 1.5 + 0.3 (EPIC) and 1.6 ± 0.2 (RGS). In 
the wgem and gdem results the relative jump in O/Fe is equal or 
lower, but still consistent with the result for the single temper- 
ature model: 1.4 + 0.3 (EPlC/wdem), 1.4 ± 0.2 (EPlC/gdem), 
1.6 ± 0.3 (RGS/wdem) and 1.5 + 0.3 (RGS/gdem). 

4.2. Abundances and SNIa/SNII/Population-lll models 

From the single-temperature (CIE) and DEM models we fit to 
the data, we obtain the abundances of the elements for which 
line emission is detected. Assuming that all the elements origi- 
nate from SNIa, SNII and PopIII stars, we can construct a sim- 
ple model to obtain the relative contribution of these objects to 
the enrichment of the ICM. 

We use several SNIa yields obtained fro m two phys- 
ically different models adapted from Jwamoto et al.l (Il999h 
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Fig. 6. RGS spectrum in the 19-21 A band featuring the O viii 
Ly-a line. The data points represent the measured spectrum, 
while the model line shows the line when the scale factor is 
frozen to ~0.3, which is the best-fit value for the iron lines. 

to fit our abundances. The W7 model describes a so-called 
slow deflagration model, while the WDD2 is calculated us- 
ing delayed-detonation (DD) models, which is the currently 
favoured Type la explosion scenario. Note that with SNll we 
mean all types of core-collapse supernovae including types lb 
and Ic. We use the SNll yields integrated over the stellar popu- 
lation calculated by Tsujimoto et al. ( 1995) and Iwamoto et alj 
For the PopIII-sta r SN yields we use two models from 
[Heser & Wooslev ('2002") with difi'erent core masses of the 
Poplll star: 65 Mq and 130 Mq. These two masses are the low- 
est and the highest core mass considered in lHeger & WooslevI 
ll2002h . 

For every element / the total number of particles A^,- is a 
linear combination of the number of atoms produced by super- 
nova type la (Y,.ia), type II (Y,- u) and PopIII stars (Y;,iii). 

A^; = aYij, + bYiji + cY,,ui, (8) 

where a, b and c are multiplication factors of SNIa, SNII, and 
PopIII stars respectively. The total number of particles for an 
element can be easily converted into a number abundance. This 
reduces to a system of three variables (a, b an c) and nine data 
points (O, Ne, Mg, Si, S, Ar, Ca, Fe and Ni). We present the 
ratio of the relative numbers of SNIa, SNII and PopIII with 
respect to the total number 

We fit the abundances we obtain from the EPIC results 
of the core region to a model consisting of SNIa, SNII and 
PopIII-star yields in order to determine their relative contribu- 
tion to the ICM. In Table Owe compare two SNIa yield mod- 
els fitted together with the SNII and PopIII yield models to 
the measured abundances in Table 0] We use the abundances 
derived from the fits with the wdem model. If we include the 
calcium abundance, none of the fits are statistically acceptable. 
However, if we ignore the calcium data point, we obtain ?lxI of 
about 1 for the WDD2 model. We also test other SNIa models 
listed in .Iwamoto et al..a999.) : W70, WDDl, WDD3, CDDl 



Table 7. Fit results for spatially resolved RGS spectra be- 
tween 8-38 A excluding CCD2. The fitted models are single- 
temperature CIE (s), wdem (w) and a Gaussian DEM model (g). 
Emission measures {Y - J n^nndV) are given in 10^^ cm""*. 
The iron abundance is fixed to 1 .0 with respect to solar abun- 
dances. 



Parameter 


Mod 


0-0.5' 


0.5-1.0' 


1.0-2.0' 


Y 


s 


3.69±0.04 


1.80+0.04 


1.46+0.03 




w 


3.60+0.04 


1.78+0.04 


1.30+0.05 







3.74±0.04 


1.88+0.05 


1.46+0.04 


kT 


S 


3.23±0.09 


4.2+0.2 


4.3+0.3 


(keV) 


5 


4.09±0.16 


5.1+0.6 


4.3+0.3 


kT 

mean 


W 


3.53±0.14 


4.7+0.4 


4.3+0.5 




w 


4.76±0.18 


6.6+0.5 


6.5+0.7 




w 


'1^+0 04 


66+0 09 

\J -i—\J .\J y 


1 06+0 18 




g 


n ^49+0 01 8 
\j .j'-TiLzn.yj .\j L 


35+0 06 

\j . J J^yJ .\J\j 


00+0 1 7 


O/Fe 




5^+0 05 


0.77+0.12 


Q6+n 1 5 






51+0 05 


77+0 1 7 


1.2+0.3 




g 




U.UO^VJ. 1 




Ne/Fe 


s 


1.10+0.12 


1.7+0.3 


0.8+0.4 




w 


0.79+0.11 


0.9+0.3 


0.00+0.10 




g 


0.75+0.12 


0.9+0.4 


0.8+0.4 


Mg/Fe 


s 


0.16+0.10 


0.00+0.10 


0.00+0.04 




w 


0.32+0.11 


0.0+0.2 


0.00+0.12 




g 


0.28+0.11 


0.00+0.11 


0.00+0.05 


Scale 


s 


0.50+0.08 


1.14+0.17 


1.3+0.3 




w 


0.45+0.07 


1.1+0.2 


2.8+0.4 




g 


0.40+0.07 


0.8+0.3 


1.3+0.2 


X^ld.o.t 


s 


1100/858 


970 / 858 


925 / 858 




w 


1000 / 855 


953 / 855 


885 / 855 




g 


1039 / 857 


971 /857 


926 / 857 



and CDD2. Their best-fit ratios were in general similar to the 
WDD2 and W7 models, but with a higher 

The relative contribution of PopIII stars is in all fits smaller 
than three times its error, thus it is not significantly detected. 
Because of the large yields per PopIII event, a small number of 
PopIII stars can in principal contribute a lot to the abundance. 
The abundance pattern of PopIII stars, however, resemble the 
patterns of SNIa and SNII depending on the core mass. The 
SNIa or SNII patterns can mostly compensate for the PopIII 
contribution, when it is left out in the fit. The does improve 
only marginally when we add PopIII star yields to a model con- 
taining only SNIa and SNII yields. The best fit using just SNIa 
and SNII models including calcium gives a;^f^/dof of 15 / 7 and 
a SNIa contribution of 0.35+0.03. In fact some models fit the 
data best by putting a negative value to the PopIII or SNII con- 
tribution, which is unphysical. The 65 Mq PopIII model pro- 
duces a sUghtly higher PopIII contribution than the one for 130 
M0. But since the mass of the expelled material is higher for 
the latter model, the number of stars needs to be higher in the 
65 M0 model to get a comparable effect. The errors on these 
bigger values are also too large to claim a significant PopIII- 
star contribution. 

In Fig. |9] we plot the abundances from the wdem model 
with the best-fit model (WDD2) listed in Table El Note 
that the WDD2 model is also the type la model favoured 
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Fig. 9. The elemental abundances derived from a wdem fit, fitted with the WDD2 SNIa-yield model. The best fit from Table|8lis 
shown in the left panel with a Poplll mass of 65 Mq. In the plot we also show the SNla, SNll and Poplll contribution separately. 
In the right panel we show a fit without a Poplll contribution. In both plots the calcium abundance is included in the fit. 



Table 8. Relative contribution of SNla, SNII and Poplll stars 
to the enrichment of the ICM. We compare two SNla models 
in llwamoto et alJ (11999.) with the data. We only show results 
from wdem fits, because we could only get an acceptable fit 
using these data. The values shown here are the fractions with 
respect to the sum of all contributions (SNla + SNll + Poplll). 



Type 




BOMq 




65Mo 




Value 




Value 


1 dof 


W7 


SNla 


0.33±0.09 




0.30±0.04 




SNII 


0.7±0.2 


All 6 


0.70±0.14 


55/6 


Poplll 


(3.4±1.4)xl0- 


-3 


(-3±7)xlO-' 




WDD2 


SNla 


0.38±0.08 




0.48±0.10 




SNII 


0.63±0.15 


12/6 


0.5±0.2 


12/6 


Poplll 


(-1.8±1.2)xlO 


-3 


(1.7±l.l)xl0 


-2 



Fits with the calcium abundance excluded: 



Type 


Mp„pin = 130Mo 




65Mo 




Value 1 dof 


Value 


;r-/dof 


W7 


SNla 


0.34±0.09 


0.32+0.05 




SNII 


0.7±0.2 29/5 


0.69±0.15 


41/5 


PopHI 


(2.9±1.3)xl0-' 


(-2±7)xl0-' 




WDD2 


SNla 


0.38±0.08 


0.50±0.10 




SNII 


0.62±0.15 4.9/5 


0.5±0.2 


4.1/5 


Poplll 


(-I.9±1.2)xl0-3 


(1.9±1.2)xl0 


-2 



bv llwamoto et al.l (Il999l) on observational grounds. The data 
points we use are taken from the third column of Table|4l From 
this plot we see that the calcium abundance might be under- 
estimated by the nucleosynthesis models and the cause of the 
high values. But the actual uncertainties in the measured 
calcium abundance might be bigger than the statistical error in- 
dicated here. This calcium overabundance is also seen in the 



analysis of 2A 0335-1-096 dWerner et al.l2006l) . confirming that 
the overabundance is probably not a statistical deviation. Every 
other abundance is consistent with the model within the error 
bars. The contribution of Poplll can be seen mainly in oxygen 
and neon, but the uncertainties are large. 

We applied the same procedure to the abundances obtained 
from RGS. The relative abundances determined from RGS are 
very well constrained, but the number of elements that we can 
measure is too small to get a reasonable fit. If we fit the abun- 
dances we obtain eiTors which are equal to the measured val- 
ues. Therefore, we do not show these fits in this paper. 

4.3. Non-thermal X-ray emission? 

In pr e vious papers about S ersic 159-03 (e.g. iKaastra et alJ 
l2003[ E onamente et alJ [2005) the larger x^- for the single- 
temperature model is attributed to a soft X-ray excess. There 
are, however, a number of solutions to fit this soft excess and to 
obtain an acceptable x]- of about 1 .0. From the spectral fits in 
Sect. 14. II we confirm that a single-temperature model is not the 
best description for the observed spectra. If we fit the data with 
a multi-temperature model, then we obtain ?l x^- of about 1.0. 
In this section we explore the possibility of the existence of a 
soft excess in Sersic 159-03 and speculate about its nature. The 
wdem and gdem models already provide acceptable fits to the 
spectra. We verified that the non-thermal component presented 
in this section does not affect the trends observed in the thermal 
analysis listed in Sect. 14. II 

One possibility is that the soft X-ray emission origin ates 
from warm thermal emission from the WHIM ( Kaastra et alJ 
l2003l) . Here we explore the possibility that inverse-Compton 
emission from CMB photons which are up-scattered to X- 
ray energies by relativistic electrons can also cause a soft ex- 
cess. The same emission mechanism is also thought to be 
the origin of the hard excess in clusters of galaxies detected 
by BeppoSAX (e.g. Fusco-Femiano et al...2005.) . A power-law 
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Fig. 10. This spectram shows that an underlying power-law component can cause both a soft and a hard excess. The figure shows 
a power-law component with F = 2, a single-temperature CIE component with a temperature of 2.5 keV and the total spectrum. 
The contribution of the power-law component is mostly noticeable below ~0.5 keV and above ~5 keV, hence causing a soft and 
hard excess with respect to the thermal component. 




Radius (arcmin) Radius (arcmin) 



Fig. 11. Left: The 0.3-10 keV integrated intensities for three model components: Cluster hot gas (□), power law (•) and soft- 
excess (*). Right: Ratio between the intensities: power law/hot cluster gas (•) and soft-excess/hot cluster gas (*). 
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Table 9. Fit results for an EPIC spectrum extracted from a 
circle with a radius of 4' and centred on the core including a 
single-temperature thermal component + power law and multi- 
temperature thermal components + power law. Fluxes are cal- 
culated over the 0.3-10 keV range and presented in 10 "* erg 
cm~^ s"' deg"^. Emission measure (Fthermai = / n^nHdV) is 
given in 10^^ cm""* and Fpow is given in 10^' ph s"' keV"' at 1 
keV. 



Parameter 


single-temp + pow 


wdein + pow 


gdein + pow 


^thermal 


16.8 ±0.3 


17.6 ±0.3 


18.9 ±0.4 


^thermal 


16.2 ±0.2 


16.9 ±0.3 


18.0 ±0.5 


kT 


2.45 ± 0.02 




2.50 ± 0.03 


kT 

"■-^ mean 




2.53 ± 0.07 








3.14 ±0.06 




a 




0.32 ± 0.04 




(Tj 






0.208 ± 0.012 


o 


0.27 ± 0.05 


0.26 ± 0.05 


0.17 ±0.03 


Ne 


1.00 ±0.11 


0.42 ±0.11 


0.11 ±0.09 


Mg 


0.18 ±0.04 


0.21 ± 0.05 


0.18 ±0.04 


Si 


0.32 ± 0.02 


0.33 ± 0.02 


0.33 ± 0.02 


S 


0.22 ± 0.02 


0.25 ± 0.02 


0.25 ± 0.02 


Ar 


0.20 ± 0.05 


0.23 ± 0.05 


0.26 ± 0.06 


Ca 


0.42 ± 0.06 


0.45 ± 0.06 


0.52 ± 0.07 


Fe 


0.433 ±0.011 


0.409± 0.011 


0.275 ± 0.009 


Ni 


0.62 ±0.11 


0.60 ±0.10 


0.54 ± 0.10 


Y 

pow 


6.6 ± 0.4 


5.3 ±0.4 


3.5 ± 0.8 


F 

^ pow 


3.5 ± 0.2 


2.8 ± 0.2 


1.6 ±0.4 


r 


2.06 ± 0.03 


2.10 ±0.04 


2.30 ± 0.15 




897/914 


853/913 


861/913 



component describing this inverse-Compton emission might be 
able to fit both this hard and soft excess. 

A single non-thermal component as illustrated in Fig. ^| 
with a r of 2.0 can explain both hard and soft X-ray excess 
with respect to a single-temperature model. At low energies 
<0.5 keV the power-law component is comparable in flux to 
the thermal component, hence causing a soft excess. Above ~5 
keV the power-law component in this example is even stronger 
than the thermal component, causing a hard X-ray excess. In 
Sersic 159-03 the power-law emission is not as strong as in this 
example, but a similar model is consistent with the data. 

In order to learn whether the data allow the soft excess to 
be thermal or non-thermal in nature, we fit the spectra from all 
annuli (see Table |3} with a single-temperature cluster compo- 
nent, a cold thermal component representing the WHIM, and 
a power-law component representing the soft and hard non- 
thermal emission. We first fit all annuli with the photon index 
and the kT of the cool component as a free parameter Because 
the photon index and temperatures were not well constrained in 
the outer annuli, we fix the photon index and cool temperature 
in the final fits of all annuli to the mean value in the core re- 
gion, which is 1 .9 for the power-law index and 0.25 keV for the 
cool thermal component. This assumption for F is valid if the 
relativistic electrons are produced by acceleration in shocks, 
which behave similarly regardless of the position in the clus- 
ter The radial profiles of the fitted components are shown in 
Fig-E] The profile of the hot gas and the soft component have 



the same shape, which suggests that they have the same ori- 
gin. The power-law component, however, has a more extended 
radial distribution. The flux is still the highest in the core of 
the cluster, but the slope of the radial profile is more shallow 
than the thermal emission. The ratio between the components 
is plotted in the right panel of Fig. ^2 This figure indeed shows 
that the ratio between the soft excess and the cluster gas is con- 
stant. This is is not what you expect if the excess is due to 
WHIM which should be picked up preferably at the outskirts 
of the cluster. On the other hand, the power-law component 
does become relatively more important in the outskirts of the 
cluster. This is to be expected for Inverse-Compton emission 
because it scales with density n instead of which is true for 
thermal emission. The different profile of the power-law emis- 
sion therefore strongly suggests that if there is an excess, it is 
non-thermal in nature. 



In the following fits we drop the cool component and fit 
the spectra with just a thermal (DEM) components and an ad- 
ditional power-law component to fit the soft-excess (Table |9}. 
The models including a power law generally result in a 
lower x^- This effect is most strongly present in the single- 
temperature fit: - 897 / 914 with a power-law component 
and x^ - 1228 / 916 without. But also for the DEM models 
there is an improvement in x^ that provides an argument for 
the presence of a non-thermal component. 

The absolute values of the abundances change when an ad- 
ditional power law is used, but the general trends do not change 
considerably. The spatial abundance distributions are also con- 
sistent with the trends observed in the single-temperature mod- 
els. In particular, fits including a power-law component show 
that the jump in 0/Fe in the centre of the cluster remains. 
However, the significance of the jump decreases compared to 
the single-temperature result when a power-law component is 
added. 



Our suggestion that a single power-law component can 
both explain a possible soft as well as a hard excess of 
Sersic 159-03, may also be relevant for several clusters for 
which non-thermal emission has been de tected by BeppoS AX 
frusc o-Femiano et al.l2000[l2004Eo05 t). However, the Coma- 
cluster is, so-far, the only cluster for which both a hard 
X-ray component and a soft excess has be en reported 
("Fusco-Femia no et al.l 12004*: 'Kaa stra et al ]L2003'). For Sersic 
159-03, it does indeed seem to be the case that a single non- 
thermal component explains both the soft and hard X-ray ex- 
cess: Extrapolating the BeppoSAX flux in the 20-80 keV 
range, (1.5±0.5)xlO"" erg s"' cm"^, to the 0.3-10 keV band, 
assuming F = 2.0, we obtain a flux of (3.8+1.3)xlO"'' erg 
s"' cm"^. This is consistent with the XMM-Newton flux for 
the putative non-thermal component reported by Kaastra et alJ 
(2003), (6.8+2.0)xl0"" erg s^' cm^-. For this calculation 
we assume that all the non-thermal emission originates from 
within a circle of 12' from the core. 
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5. Discussion 

5.1. Temperature structure 

The parameters of the DEM distributions which we obtain from 
multi-temperature fitting, are well determined for this cluster. 
However, we cannot discriminate between different shapes of 
DEM distributions and the presence of a non-thermal compo- 
nent. The DEM parameters a (wdem) and (Tj (gdem) values 
show a slight increase to the outer parts of the cluster. This can 
be explained by the steep temperature gradients we see both in 
the cooling core and in the outer parts of the cluster The val- 
ues we der ive for a are higher than the value of 0.20 + 0.05 
derived by iKaastra et alJ ( l2004l) based on a shorter exposure. 
However, their background subtraction method and handling 
of spectral excess is different, which can lead to systematic dif- 
ferences. Also projection effects can hav e an influence on the 
broadness of the DEM distribution, but as'Kaa stra et alJ ll2004!) 
and Werner e t al. ( 2006 ) show, multi-temperature models are 
also needed for fitting deprojected spectra. 

5.2. Abundance distribution and enriciiment by 
supernova types la/ll and Population-Ill stars 

From the EPIC and RGS spectra we obtain radial abundance 
profiles for the most abundant metals. The EPIC and RGS ra- 
dial profiles for oxygen show a jump in the centre of the cluster. 
The key data point is the point in the 0-0.5' bin. If we compare 
the O/Fe ratio of the central bin with the neighbouring bins, 
its value is significantly lower This means that either the iron 
abundance is relatively high or oxygen is low in the centre of 
the cluster. 

Recently, using cluster simulations, ISchindler et al ] il2005l) 
found that ram-pressure stripping acting on in-falling galaxies 
can result in centrally peaked abundance profiles, while early 
galactic winds produce an extended distribution of the oxygen 
abundance. The spatial distribution of oxygen, which is difficult 
to determine, is the key to understand the evolution history of 
the cluster 

Contrary to the iron and oxygen lines in Sersic 159-03, the 
widths of these lines are the same in the cluster 2A 0335-1-096 
jWerner et al. 2006 ). A possible explanation is that the clus- 
ter 2A 0335-H096 shows a much stronger temperature drop in 
the centre of the cluster. Since spectral lines are stronger when 
the temperature is lower, the RGS spectrum is dominated by 
line emission from the cool core of the cluster In this case, the 
line profiles of iron and oxygen should follow the temperature 
structure more than the abundance distribution of the elements. 
In Sersic 159-03 the temperature profile within the RGS extrac- 
tion region is relatively flat. The temperature is about 2.5 keV 
with a spread of 10%. Because a strong thermal gradient is ab- 
sent in the core of this cluster, the line profile should follow 
the abundance distribution of the elements. The observed line 
profiles in RGS therefore strongly suggest that in Sersic 159- 
03 oxygen emission is lower in the centre of the line profile 
compared to iron, which is consistent with the radial profiles 
extracted from EPIC and RGS. 



From several abundance yield simulations for Supernova 
type la and II we know that oxygen dominates the yield of type 
II supernovae, while iron originates mostly from type la super- 
novae. If the iron abundance is more peaked in the centre of 
the cluster, we can formulate a scenario of the possible enrich- 
ment history of the ICM in Sersic 159-03. In the early universe 
the Inter-Galactic Medium (IGM) might have been enriched by 
PopIII-star explosions. The star bursts in the young galaxies 
that form after PopIII enrichment produce a lot of SNII which 
produce mainly oxygen. Because clusters are still in an early 
stage of development, the oxygen becomes well mixed through 
the IGM. Then, 1 billion years after the 'Big Bang', also SNIa 
explosions start enriching the local ISM in galaxies continu- 
ously. If a galaxy containing both SNIa and SNII products falls 
into a cluster, the ISM is stripped off by ram-pressure strip- 
ping, preferably in the denser core of the cluster, hence lead- 
ing to a more centrally peaked abundance dis tribution of all 
SNIa and SNII products ( Schindl er et alJl2005l) . Because iron 
is also produced by SNIa explosions in the central cD galaxy, 

this could enhance the i ron abundance in the centre of the clus- 
j ■ I 

ter ffamura et al. 2001). The drop in the centre of the observed 
O/Fe fits within this scenario. 

In order to check whether the absolute abundance of iron 
can be consistent with the interpretation we just formulated, 
we calculate the total iron mass in the cluster We find about 4 
X 10*^ Mq within a radius of 300 kpc of the core. If all this iron 
originates from SNIa explosions, we would need at least 5 x 
10'^ explosions in the cluster lifetime, which is of the order of 
10"^ years. Since there are typically 10^ galaxies, there should 
be about 1 SNIa explosion per 200 years per galaxy, which is 
consistent with the rate of SNIa in our own galaxy. 

We then fit the abundances from E PIC to a model of the 
yields of S NIa, SNII and PopIII stars ("iwamo to etaDll999t 
iHeger & Wooslev 2002) . Using a linear combination of the 
yields we can estimate the expected abundances and fit them 
to the data. From the fit we obtain the relative contributions of 
the three types of metal-enriching sources. 

We have to be careful, though, when we interpret the fit- 
ted supernova ratios. The mechanisms by which the elements 
are ejected into the ICM are not yet fully understood and are 
certainly not included in our simplified model. A small part 
of the supernova products can also be locked up in low-mass 
stars, so the abundance distribution we observe in the ICM 
is not entirely representative for the total abundance distribu- 
tion in the cluster. A more detailed description of the prob- 
lems when one uses this sim plified model is described in 
iMatteucci & Chiappini ll2005l) . The derived supernova ratios 
should be interpreted as the number fraction of supernovae that 
would be needed to enrich the ICM, not the actual number of 
supernovae during the history of the cluster. 

Our analysis shows that PopIII stars are not required to 
fit the data. The contributions from this early generation of 
stars are probably too small to be detectable in current clus- 
ter abundance patterns. We cannot confirm the results from a 
sample of clusters observed with ASCA that suggest that low 
mass PopIII stars are nec essary to explain the o bserved sili- 
con and sulfur abundance (Baumgartner et al. 2005). Our anal- 
ysis, which includes more elements with about the same error 
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bars on the abundances as in lBaumgartner et alJ ( l2005t) . shows 
that a PopIII-star contribution is not detected in Sersic 159- 
03 . Considering the s ingle-temperature spectral models used 
bv lBaumgartner etaOIZQOS), we think that a temperature bias 
in the iron abundance makes their Si/Fe and S/Fe ratios unre- 
liable. We also ca nnot confirm the u nderabundance of argon 
and calcium in the iBaumgartner eTal j (^2005) sample. We ac- 
tually measure an overabundance of calcium with respect to 
the nucleosyn thesis models, which is in line with the result by 
IWerner et all (t2006 ) in the cluster 2A 0335 -h096. This might 
be an indication that the SNll and SNla models have difficul- 
ties predicting the right yield for this element, but the actual 
uncertainties in our measured value might also be bigger than 
the quoted statistic al error. When we c ompare the supernova 
fractions found bv IWerner etaP ilOO ^) in 2A 0335+096, we 
see that the SNla contribution in Sersic 159-03 tends to be 5- 
25% higher 

5.3. Soft excess and non-thermal X-ray emission 

An alternative solution to obtain an acceptable fit for the spectra 
of Sersic 159-03 is to add a soft-excess component to a single- 
temperature model. We explore the possibility that the devia- 
tions from the single-temperature model are due to a different 
component, Uke the WHIM or inverse-Compton emission. 

A non-thermal component can explain both a soft and 
hard excesses for both Sersic 159-03 and Coma, is consistent 
with inverse-Compton emission of CMB photons on relativis- 
tic electrons (iHwanglll997h . The alternative explanation, non- 
thermal bremsstrahlung emission from a non-thermal tail to the 
Maxwellian energy distribution of the electrons, would only re- 
sult in a hard X-ray excess. 

There are several possible sources that can produce a pop- 
ulation of relativistic electrons in a cluster that would explain 
the observed power-law emission. The acceleration of electrons 
can occur, f or example, in active galactic nuclei, supernovae 
and pulsars dSarazin & LieuiJ99 8). But also plas ma waves or 
shocks in the ICM can accelerate electrons (see lEnBlin et al.l 
1 1999 for a discussion). Using our results from Sersic 159-03, 
we propose a new scenario where electrons are accelerated by 
the shock waves associated with in-falling galaxies. 

We already discussed the possibility that the centrally 
peaked iron distribution might be due to an eflicient ram- 
pressure stripping of in-falling galaxies in the dense core of the 
cluster Some fraction of these galaxies being stripped would 
have moderate supersonic velocities. The galaxies will drive 
bow shocks of Mach numbers up to ~ 3 in the intra-cluster gas 
that can accelerate electrons to relativistic energies. To produce 
photons of ~ 10 keV due to inverse-Compton up-scattering of 
CMB photons, electrons with a Lorentz factor y ~ 7 x 10^ are 
required. The life-time of such an electron in the cluster core 
is about 2.3xlO'V ' years (assuming that the magnetic field 
in the core is below 2/iG and the energy densities of optical 
and IR radiation field in the core are below 0.1 eV cm"^). It 
is worth noting here that the central galaxy in Sersi c 159-03 is 
one of the brightest cD galaxies detected in IR (.Hansen et alJ 



I2OOOI) . but the energy density of the IR photon field produced 
by the galaxy is below 0.1 eV cm"^. 

The electron acceleration time at a shock of velocity > 1000 
km s"' in the cluster core is typically much shorter than 3x10*^ 
years. Thus, to provide a steady non-thermal X-ray emission 
the stripping galaxies must cross the cluster core at a frequency 
V of about 1 in 3x10^ years. The mechanical power, Lyn, dissi- 
pated by the galactic shocks of velocity Vg and average radius 
Rg in the core region of diameter and average number den- 
sity He can be written as: 




For the parameters assumed above, galaxies with velocity Vg ~ 
2000 km s"' would provide L^in ~ 10'*'* erg s"'. A few per- 
cent of the power would be enough to maintain the relativistic 
electrons required to provide the 10 keV regime non-thermal 
emission. Note that the lifetime of the electrons with Lorentz 
factor below 100 would be about the Hubble time. This implies 
that the spectrum of electrons accelerated by multiple succes- 
sive shocks of Mach number about 3 will have power-law dis- 
tribution of index about 3 at 7 ~ 10^, and a flat distribution 
at 7 < 100, where the Compton losses are not important. The 
electron distribution index 3 coiTesponds to the photon index 
of the photon emission of 2, close to that observed. 

6. Conclusions 

We have analysed high-resolution X-ray spectra of the clus- 
ter of galaxies Sersic 159-03 obtained with XMM-Newton and 
conclude that: 

1 . For the first time we accurately determine the abundances 
in a cluster of galaxies which fit to a linear combination of 
current supernova yield models. 

2. From the line width in RGS and the radial profiles from 
EPIC/RGS we find a jump in the O/Fe ratio around a ra- 
dius of 0.5' from the cluster centre. The O/Fe ratio in 
the centre of the cluster is lower compared to its immedi- 
ate suiToundings. A co mbination of ram-pressure stripping 
dSchindler et al J2005l) and enrichment by SNla in the cen- 
tral elliptical galaxy dTamura et alJbOoH) can explain the 
observed profile. 

3. We do not detect a contribution of PopIII stars in Sersic 
159-03. This result is not in line with the claim by 
IBaumgartner et all (l2005h that Poplll stars are necessary to 
explain the abundances measured in an ASCA sample of 
clusters. The number contribution of SNla with respect to 
the total number of supernovae we find, based on the mea- 
sured abundances, is about 25-50%. 

4. The spectra can also be fitted with an additional non- 
thermal component with a power-law index of about 2. 1 on 
top of the thermal emission. The of these fits is compa- 
rable to the using only multi-thermal models. If a non- 
thermal component is present in Sersic 159-03, it can be 
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explained by Inverse-Compton scattering of CMB photons 
and relativistic electrons. 
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